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H306peTeHHe OTHOCHTCH K KOHTpOJIb- 
HO-M3MepHTeJIbH0H TeXHHKG H M03KGT 6bITb 

HcnonbsoBaHO win HSMepeHHH pasMepa 
KaneJib ocaAKOB c noMombio ynbTpasBy- 
. KOBblX Kone6aHHH. 

UeJIbK) H306peTeHHH HBJIH6TCH nOBbl- 

meHHe to^hocth HSMepeHHH. 

Ha MepTesce npeACTaBJieHo ycTpoacT-r 
bo ajih peajiHsanHH cnoco6a HSMepeHHH 
pa3Mepa Kanejib oca^KOB. 

Yctpohctbo conepscHT renepaTop 1 » 
coeflHHeHHHH c aKycTH^ecKHM Hsny^a- 
TejieM 2, Kopnyc 3, b kotopom noMeme- 



HSMepeHHH. CnocoS HSMepeHHH pasMepa 
nanenb ocajjKOB ocymecTBJjHeTCH cneAy"" 
wmHM o6pasoM. B Hccjie/xyeMOH cpeAe 
BO36yroai0T ynpyrae KOJie6aHHfl. IIphhh- 
MawT aicycTOTecKHH cKrHan, npomewnHH 
3a,aaHHoe paccTOHHHe b HCcneAyeMon 
cpefle. HsMepHioT aMruiHTyAy npHHHToro 
cnrnana. IIocKonbKy npH HaxomAeHHH 
KaruiH b 3one pacnpocTpaHeHHH ynpyrux 
bojih tiacTb sHeprHH oTpaxaeTCH ot *iac- 

THUH oSpaTHO, TO H aMtUIHTyAcl npHHK- 

Toro aKycTH^ecKoro cnrHajia 6yneT • 
yMeHbmaTbCH, T.e. sa BpeMH npeSbiBa- 
hhh Kama* b To^xe axy CTH^iec khx bojih 
6yAeT npoHcxoAKTb MOfly jihahh iiphhhto- 
ro aKycTH^tecKoro cHrHajia. MeM 6ojibme 
tiacTHixa, TeM c 6ojibmeH cKopocTbio oHa 
nawaeT h TeM MeHbme 6yner BpeMH mo~ 
jxynxiXKH aMnnHTyAbi. H3MepHH)T BpeMH 
MOAyjiHUHH aMiumTyAM npHHHToro aicyc- 
TtmecKoro cHraajia h no H3MepeHH0My 
BpeMeHH onpeAeJiHioT pasMep KaneJib 

OCaflKOB. 1 HJl. 
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Hbi aKycTHtiecKHH H3jry^aTejib 2 h aicyc- 
THMecKHH' npHeMHHK 4. K Bbixofly npneM-' 
HHKa 4 nocjieAOBaTejibHo noAKJiroueHbi 
ycHjiHTejib 5 f AeTeKTop 6 h ceneKTop 
7 AJiHTeJibHocTH. B BepxHeft ^acra Kop- 
nyca 3 HMeeTCH oTBepcTHe 8. Buxoah 
: cejieKTopa 7 A-nHTenbHocTH coeAHHenbi 

.C COOTBeTCTByKmHMH BXOAaMH 6jioKar .9 
CUeTHHKOB . 

CnocoS HSMepeHHH pasMepa xanejib 
ocaAKOB ocymecTBJiHWT cneAywmHM o6pa- 
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B K.opnyce 3 co3flae*rcff c noMombw 
axycTHwecxoro H3Jiy*iaTejiH 2 axycTH- 
qecKoe noJie. AxycTH**ecKHH H3JiyHaTenb. 

2 B036y»maeTCH ot reHepaTopa 1 3nex- 
TpH^ecKHx CHrnanoB, CosAaBaeMoe BHyT- 
PH Kopnyca 3 axycTHMecxoe nojie npHHH- 
MaeTCH aKycTH^ecKHM npHeMHHKOM 4 , 
ITpeo6pa30BaHHbift CHrHan 3aTeM nocTy- 
naeT Ha ycwiHTejib 5, AetexTop 6 h 

CeneKTOP 7 AJJHTeJIbHOCTH. IIpH OTCyTCT"" 

bhh *iacTHU. ocaAKOB b Kopiiyce 3 cnr- 
Han, npHHHMaeMH# axycratiecxHM npHeM- 
HHKOM 4, He HSMeHHeTCH no aMnjiHTyae. 
HacTHUbi ocatfKOB, nonaAawmne b KOpnyc 

3 ^epe3 oTBepcrae 8, oTpascaioT *tacTb 
aKycTH^ecKott 3HeprnH t Bp36y?KAaeMOH 
aKycTHMecKHM H3JiynaTejieM 2. IIpH stom 
yMeHbinaeTCH aMruiHTyAa axycTH*tecxoro 
cnrHajia, npHHHMaeMoro npneMHHKOM 4,. 
J3to H3MeHeHHe aMruiHTyAbi HaojnoAaeTCH 
b Te^eHHe Bcero BpeMeHH npedbisaHHH 
KaiuiH b Kopnyce 3 npH ee najjeHHH, 
T.e. npH npoxoH«eHHH ^acTHi^i ocaAXOB 
^epe3 Kdpnyc 3 npcmcxoAHT MoayjiHUHH 
aMiuiHTyflfci npHHHMaeMoro axycTH^eckoro 
cHrHana* CKopocTb nafleHHH xaiuiH oa~ 
H03Ha^Ho cBH3aHa c ee pa3MepaMH. 3to 
nosBOJineT no H3MepeHHofi ckopocth na- 
AeHHH KaiuiH oueHHTb ee pa'3Mepu« IIpH 

3TOM B SaBHCHMOCTH OT CKOPOCTH H3Me~* 

HjreTCH BpeMfl npeobreaHHH tiacTHUbi b 
Kopnyce 3. TaKHM o6pa30M, HH$opMaijHH 
o pa3Mepax xamiH nony^aeTC* He no 
napaMeTpaM npHHHToro axycTH^tecxoro 



cHrHajia, a no fcaxTy HariHMHH MpnyjiH- 
Uhh 3Toro CHrHajia b TeneHH onpeAe- 
jieHHoro BpeMeHH, mto AaeT bosmokhoctb 

5 nOBblCHTb TOMHOCTb HSMepeHHft. H3MepH- 

eTCH BpeMH, b Te^eHHe KOToporo H3Me- A 
HHeTCH aMnjiHTyAa npHHHToro axycTH- 
tiecKoro cHrHana, T,e« npoHcxoAHT ero 
MOAynHUHH . Jinx 3Toro ceneKTop 7 ajih- 
10 TeJibHocTH ocymecTBJiaeT cejiexuHK) npH>- 

« HHTbK HMnyJIbCOB, HanpaBJIHH HX B Sa- 
BHCHMOCTH OT AHHTeJIbHOCTH Ha OAHH 
H3 CHeTHHKOB 6nOKa 9 C^ieTTOXOB. Ko- 

nHMecTBO CMeTMHKOB o npe Ae nne tc h He- 
15 odxoAHMOH cTeneHbw tomhocth HSMepe- 

HHft. KaKAblH H3 CMeTMHKOB 4>HXCHpyeT 

HajiHMHe MacTHObi c pasMepaMH, cootbct- 
cTBywmHMH onpeAeneHHOMy AHanasoHy, 

2o<J>ppMyjia H3o6peTeHHH 

Cnoco6 H3MepeHHH pasMepa xanejib 
ocaAKOB, ocHosaHHbifi Ha HanpaBneHHH 
Ha ocaAKH nyqxa axycTHMecKHX xone6a- 

25 hhh, npneMe sthx xoJie6aHHft, HSMepe- 
hhh aMnjiHTyAbi npHHHToro cHraana h 
onpeAeneHHH pa3Mepa xanejib ocaAKOB 
c y^eTOM H3MepeHHOH aMruiHTyAM, o t - 
ji h m a H) m h ft <T h TeM, MTO f; c ue- 

30 Jibio noBbnneHHH tohhocth HSMepeHHH, 
nponycxaioT xaruiH ocaAxoB tiepes nywx 
axycTHMecKHX Kone6aHHtt f HSMepnioT Bpe- 
mh H3MeHeHHH npHHHToro CHrkana h on- 
peAenHioT pasMep xanenb ocsakob no hs- 
.MepeHHOMy BpeMeHH • 

*35 .! • 
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A Stochastic Raindrop Time Distribution Model 

J. LAVERGNAT AND P. GOLE 
Centre d 'Etude des Environnements Terrestre et Planetaires, Velizy, France 
(Manuscript received 25 February 1997, in final form 5 December 1997) 

ABSTRACT 

A disdrometer simultaneously measuring time of arrival and size of raindrops was set up in the Paris, France, 
area. Data collected over a period of 25 months (May 1992 to May 1994) are presented and analyzed to derive 
a long-term temporal model governed by a renewal process whose survival law is a Bi-Pareto law of the third 
kind. The model t hus found allows nearly nine orders of magnitud e of the time intervals be tween raindrops to 
be mathematically represented at the sam e time using only six parameters7 ~~ ~ 

The analysis presented here does not consul lalrtlail intensity and the nature of rain (convective, stratiform, 
etc.) as classification parameters. This approach, which may at first sight seem objectionable, is justified by the 
quality of the statistical inferences that can be made from the model. Two such applications are described — 
namely, the prediction of the total fallen-water height and the conversion between various rain gauge integration 
times, which are often necessary for telecommunications purposes (for which only limited models are currently 
available). Since this kind of temporal data is rare, a comparison is also made with published data having the 
finest possible temporal resolution from the point of view of the fractal properties of rain, namely, its fractal 
dimension. A fairly good agreement was found with these other results and at the same time leads to a deeper 
insight into the fractal nature of rain. 

This model provides a very satisfactory statistical representation of rain but does not intend to provide a 
physical interpretation of the observed temporal behavior of rain, which remains to be done. 



1. Introduction 

Rainfall is a very complex, naturally occurring phe- 
nomenon, many of the characteristics of which strongly 
influence in various ways a number of areas as different 
as hydrology (e.g., erosion and infiltration), meteorol- 
ogy, the earth *s water cycle, remote sensing (e.g., ob- 
servations from satellite platforms and radar observa- 
tions), radio communications (e.g., rain attenuation), 
etc. 

Rain is generally measured by means of rain gauges, 
which provide the total height of water fallen during a 
given time period, such as daily or hourly heights. Al- 
though this information is appropriate for climatologic 
or general meteorological purposes, its time resolution 
becomes insufficient when one tries to correlate rain 
with fast-changing phenomena, such as radar returns in 
rain or rain-induced attenuation of microwaves in tele- 
communications links. While rain gauges provide 
"flux" or intensity information — that is, an amount of 
water (height in millimeters) per time unit (hour) — an- 
other instmment^the^di ^rnrneter, measures the times 
of arr i v^L-an^-stre ^nf^ raindrops. 

The idea that rain cannot be modeled appropriately 
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by resorting only to the first group of parameters 
(fluxes), excluding the second one (sizes and times of 
arrival), can first be attributed to Smith (1993). There 
are many investigations of rain-gauge data, in which 
solutions are devised to overcome the problem of rain- 
gauge integration time (e.g., Waymire and Gupta 1981). 
For example, Olsson et al. (1993) devised a tipping- 
bucket rain gauge with a small bucket volume to obtain 
a resolution as fine as 0.035 mm min~ l . 

Most studies devoted to rain use various classifier 
parameters (Laws and Parsons 1943; Hosking and Stow 
1987). Some of these classifiers are related to the nature 
of rain, such as the notions of convective or stratiform 
precipitation. Others are related to precipitation inten- 
sity, rain rate being the most widely used. As an ex- 
ample, Feingold et al. (1986) classified raindrop size 
distributions according to rain rate. 

This paper aims at analyzing the fine temporal struc- 
ture of rain, rather than its intensity aspect, by using 
the data provided by a disdrometer that was part of a 
measurement campaign aimed at investigating the ef- 
fects of the atmosphere on millimeter wave propagation 
(Gloaguen and Lavergnat 1995; Lavergnat et al. 1988). 
This device, in addition to providing 1-min size and 
velocity dual distributions — as was its first purpose — 
was improved so as to simultaneously record the size 
and time of arrival of raindrops falling through its col- 
lecting surface. The results presented below were ob- 
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parallelepipedic light beam (250x40x10 mm) 

Fig. 1 . Schematic layout and dimensions of the optical 
disdrometer. 



tained during a total measurement period of nearly 2 yr 
for the distributions, and 4 months for the time-of-ar- 
rival data. 

The approach used in this study was to model the 
long-term behavior of rain by analyzing the whole set 
of data without attempting to classify it according to 
event types or rain intensity or any other classifier. Other 
studies that take into account preclassifications are under 
way and will not be described here. 

Nevertheless, other improvements seem to be indis- 
pensable and will be reviewed as a conclusion. The 
obtained results are encouraging and shed enough new 
light on rain structure that they were deemed useful to 
publish now. 

After a brief description of the disdrometer used in 
the experiment, the obtained measurements are pre- 
sented and a tentative interpretation is suggested in 
terms of a point process. Some interesting applications 
of the model are then derived. In the last section the 
same data are analyzed in terms of fractals and com- 
pared with already known results. 

2. The experimental setup 

The disdrometer used in the present experiment em- 
ploys an infrared beam (Hauser 1 984). An optical device 
transmits a parallel light beam (Fig. 1). On the receiving 
side, a photodiode delivers an electric signal that is pro- 
portional to the power of the received beam. 

When a raindrop crosses the beam, it causes a change 
in the diode output voltage. Each drop generates a neg- 
ative pubejikejyay_ej]p^^ the magnitude and duration 
of which are related to the size and fall velocity of the 
drop, respectively. The device can distinguish two con- 
secutive drops from each other ifj he falling edge of the 
pulse gener ated by the first on e can be Tistin^uliJrTed 
from^he_rising-edge,corresponding to the second one. 
For this to happen, the capture area has been designed 
to be sufficiently small (100 cm 2 ) to make coincidence 
of several drops in the beam negligible. Whenever co- 
incidence occurs, this leads to unrealistic velocity-di- 
ameter combinations that are easily removed from the 
obtained data. The output of the photosensor was dig- 
itized and processed by a microcomputer. A calibration 



procedure consisting of moving a rod of a size corre- 
sponding to a 2-mm drop across the beam was auto- 
matically performed once a day in order to check a 
possible drift in the device's response. Two sets of data 
were recorded: one was aimed at assessing the fine tem- 
poral or short-term structure of rain and a second one, 
with a coarser resolution, was originally intended as 
auxiliary data to retrieve rain parameters influencing 
microwave propagation in the frame of the Olympus 
propagation program (Hughes 1993). 

The first set of data was obtained by processing the 
digitized samples to provide the wanted parameters: 
1 time of ar rival t (resolution: 1 ms) and diameter D (res- 
olution: Tff^^Tnrrrr from 0 to 4.5 mm). In this dataset, 
fall velocities V were available but were used only for 
the consistency check mentioned in the previous para- 
graph. They do not intervene in the development of the 
model. This first dataset covered a continuous time span 
of 4 months, from 22 October 1993 to 20 February 1994. 
A 1-week period has also been recorded from 19 to 26 
July 1993 for a first assessment and is shown in Figs. 
2a and 2b. Note that the beam thickness (1 cm) makes 
the time separation between raindrops unreliable for 
very small raindrops (<0.5 mm) since at a 1 m s _1 fall 
velocity a raindrop resides only 10 ms in the beam. To 
overcome this problem, a new apparatus is under de- 
velopment, as the small size domain will be shown to 
be of importance here. The hardware section was de- 
veloped by the Centre d'Etudes des Environnements 
Terrestre et Planetaires in Velizy, France, while the ac- 
quisition software was developed by the Laboratoire de 
d'Etude des Transferts en Hydrologie et dans 
TEnvironnement, in Grenoble, France. 

Thus, a time series of triplets (t, D, V), characterizing 
each of the raindrops that have crossed the 100-cm 2 light 
beam, makes up the first or short-term dataset and will 
be used to measure interdrop time intervals with a 1- 
ms resolution. Though it may be objected that the time 
and space aspects of the measurement are mixed, the 
former can be considered as prevailing because of the 
small interception sectional area. 

The data were then checked for consistency. Some 
peaks occur in the background noise of the measured 
signal, especially during norirainy, fine weather periods. 
Due to their strong correlation with daily temperature 
variations, it was assumed that some of these spurious 
signals were caused by heat-induced turbulence or stray 
light. These peaks may be mistaken with pulses caused 
by actual raindrops, but on diameter-velocity (D-V) 
graphs they are easily detected as "false alarms," that 
is, drops with small diameters and very large velocities. 
To make sure these false alarms were removed, a thresh- 
old of 0.3 mm was applied to the obtained diameter 
data. It is shown in appendix A that this truncation has 
no impact on the model. This high time resolution da- 
taset spanned a total time period of 14 months, with 
several interruptions ranging in length from 1 week to 
3 months. The longest available continuous acquisition 
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period was 4 months long and is the one used in the 
derivation of the model described here. 

The second or long-term dataset was generated by 
software designed to give 1-min D-V dual histograms, 
recorded as two-dimensional arrays, and was used to 
count the number of nonrainy time intervals with a 1- 
min resolution. It spanned a 2-yr period from May 1 992 
to May 1994 (25 months). To remove the influence of 
the above-mentioned spurious signals in the analysis of 
this second dataset, only those arrays with drop counts 
greater than 60 were kept for further analysis. This se- 
lection criterion resulted in a rejection percentage com- 
parable to the 0.3-mm rejection threshold used for the 
first dataset. This procedure is justified by the fact that 
the time intervals that are looked for with this second 
dataset — namely, those greater than 10 min or so — cor- 
respond to the rain recovery time. 

3. Observations 

A first measurement campaign was carried out in 
Gometz-la-Ville, France, near Paris during a stormy pe- 
riod of 10 days (6-15 July 1993), which is not part of 
the first dataset mentioned above. Figure 2a shows a 
histogram (circles) of the decimal logarithm of time 
separations between raindrops, or interdrop intervals, 
designated r, as the random variable. This logarithmic 
random variable log 10 (r) was counted in the range from 
1 to 10 within bins of width 0.1. Note that the logarithm 
was used as the random variable due to the very wide 
range of orders of magnitude of the data (from 1 ms to 
more than 1 day). Another important point is that all of 
the interdrop time histograms presented here correspond 
to the specific measuring cross-sectional area of the dis- 
drometer used in this study (100 cm 2 ), which therefore 
conditions the obtained histograms. This difficulty is 
inherent to any local device that is aimed at measuring 
rain and always has a finite size. Although it is beyond 



the scope of this paper to address this problem more 
deeply, simulations not presented here have shown that 
the size of the capture area mainly affects the shape of 
the logj 0 (f) histograms in the lower portion of their 
range (1-4), as longer time intervals correspond to time 
separations between individual rain events, which can 
be considered independent of the measuring cross sec- 
tion. 

Figure 2b shows, as a solid line, the histogram, des- 
ignated n D , of drop diameters obtained during the same 
10-day period, from 0 to 5 mm with bin widths of 0.5 
mm. For comparison purposes with conventional drop 
size distributions (e.g., Marshall and Palmer 1948), Fig. 
2b also shows as circles the corresponding number den- 
sity of raindrops in m~ 3 mm" 1 , N D , obtained from the 
histogram, n D> through the following conversion: 

n D 

N ° ~ V(D)ATkD' 

where V{D) is the terminal velocity of a drop of diameter 
D, A is the cross-sectional area of the device, T is the 
total time during which the distribution is measured, 
and M) is the histogram bin widtfcS For V(D), the fall 
velocity relation suggested by Atlas et al. (1977) has 
been used. 

One is first surprised when looking at Fig. 2a that the 
distribution decreases rather steadily over more than 
seven decades, except for a stronger fluctuation and a 
slight change in slope after the 3-4-min histogram bin 
[logio( T ) = 3-4] due to the relatively smaller number 
of data samples above this time range. The second sur- 
prise arises from the size histogram that, in addition to 
presenting a very high regularity, shows a steep increase 
in the region of small raindrops. This observed regu- 
larity must be attributed to the 1-week integration time 
used in this particular analysis and was not specific to 
these data. 
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Fig. 3. (a) Distributions of interdrop time intervals and (b) sizes. High-resolution data from 22 October 1993 to 24 February 1994. 



Figure 3 a displays a histogram of the log 10 (r) variable 
(circles) obtained (with the same range and bin width 
as Fig. 2) for the 4-month period of uninterrupted mea- 
surements provided by the first high-resolution dataset. 
It has a shape roughly similar to that shown in Fig. 2 
[i.e., a maximum at log 10 (r) = 2, the same decreasing 
slope for 4 < log 10 (r) < 6]. It also has an even greater 
regularity than that depicted in Fig. 2, and the largest 
nonrainy interval reaches a value of a little more than 
2 days. In this histogram, the change in slope again 
clearly appears near the 5-min histogram bin. 

Figure 3b shows, as a solid line, the drop size his- 
togram obtained for the same 4-month high-resolution 
dataset. This figure also shows as circles, as in Fig. 2b, 
the size distribution of raindrops in m~ 3 mm -1 , with a 
corresponding scale on the right of the graph. For sim- 
plicity, the models developed in the following sections 
were derived from the size histograms in counts per bin 
and not from the size distributions in ntr 3 mm" 1 . This 
distribution also has a more regular shape than the 1- 
week drop size data due to the increased number of 
samples. Although less pronounced than in the 1-week 
data, it still has an increasing slope near the lower limit 
of the analyzed diameter range (0.3-5 mm). 

The possible causes of instrumental errors (e.g., spu- 
rious peaks in the recorded signal and splashing of 
drops) or software errors (for instance, counting two 
drops when there is only one, close to the background 
noise) that might have influenced the obtained drop size 
histogram shape were ruled out to explain the "abnor- 
mally" high proportion of small drops at the point where 
the apparatus reaches its sensitivity limits. The effect 
of spurious signals has been addressed by rejecting 
drops below 0.3 mm. The design of the disdrometer as 
a cylinder precludes any splashing of drops (Hauser et 
al. 1984). Multiple counting of drops leads to two con- 
secutive drops in the measured time series having the 



same time stamp, so that both of these can be rejected, 
which has been performed. 

4. Interpretation in terms of renewal processes 

The histograms shown in Figs. 2 and 3 inherently 
suggest that the arrival of raindrops may be governed 
by a point renewal process. Their regularity and repro- 
ducibility strongly favor such a process. In addition, the 
correlation coefficient between two consecutive arrivals 
computed over the time period of Fig. 2 was found to 
be 0.08, a value that may be considered characteristic 
of renewal processes that are fundamentally not corre- 
lated. 

Many authors (Waymire et al. 1981) who have in- 
vestigated this subject had rejected the hypothesis that 
rain behaved according to a renewal process because 
this was difficult to combine with the clustering behav- 
ior also obvious in rain data. Complex processes — gen- 
erally obtained by a combination of several simpler pro- 
cesses — were then suggested, such as Neyman-Scott or 
Davis (Waymire et al. 1981; Smith and Karr 1983). 

Recently, Smith (1993) proposed a marked point pro- 
cess that involves a time-varying Poisson process 
"marked" by a lognormal drop size distribution. We 
shall also assume here that raindrop arrivals follow a 
renewal process and shall try to develop an empirical 
model to represent such a process. Moreover, the clus- 
tering behavior depends critically on the survival law 
of the process. Before going into details, the employed 
methodology will be reviewed briefly. 

It should be pointed out that a renewal process may 
be determined either through the counting properties or 
the survival law. The former are simpler to obtain ex- 
perimentally, but it is generally difficult to derive the 
underlying survival law. In the present experiment, we 
had the possibility of having many observations of the 
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survival law. This is why we chose to determine this 
law first. 

When trying to determine the form of the renewal 
process in question, it may seem that a goodness-of-fit 
criterion would not be sufficient to correctly validate 
the model that is sought; however, it is a good indicator 
and should be understood as a first step in the model 
determination. The applications described in later sec- 
tions were developed in part to check the validity of 
some of the consequences of the model derived in the 
following section. Furthermore, although this antici- 
pates some of the results described below, the fact that 
this model allows the temporal distribution and drop 
size spectra to be intimately related to each other seems 
to be a very good indication of the interest of such a 
model. 

a. The renewal process 

As a first step, it will be assumed that the process is 
simple and that the first ar rival time in the r nprlplpri-Hmp 
se ries is 0. If t, = t i+[ — L is the time interval se parating 
t wo consecutive raindrops, o r interdrop interval, tne as- 
sumption of a renewal process states that the t/s are all 
independent and identically distributed. Their common 
distribution law, also known as a survival law, is either 
characterized by a cumulative distribution function F(r) 
such that 

F(t) = Prob{T, < t} 
or by a probability density function if it exists: 

dF(r) 



fir) = 



dT 



Noting that a renewal process with an exponential 
survival law is a Poisson process, the histograms found 
from the measured data described in the previous section 
show that the process we are looking for is not Poisson 
distributed, as its slope is not exponential for large r's, 
but one that varies, as a first approximation, as a de- 
creasing power of r. 

Among known probability density functions, Pareto 's 
group exhibits a power behavior for large -r's and is 
therefore a good candidate. A closer look at the data in 
Figs. 2 and 3 shows that the histograms indeed fall off 
along a straight line after having reached their peak 
around the 2.5-ms abscissa and that the change in slope, 
which may be observed around the 5-min abscissa, 
seems to indicate that a more refined distribution func- 
tion is needed. We thus attempted to derive a more 
elaborate model, though one still based upon the Pareto 
distribution. 

As suggested in the previous section, the obtained 
histograms simultaneously represent widely different 
timescales. The approach chosen to refine the initially 
selected Pareto function, suggested by the first short- 
term analyses presented in the observations section, was 
to consider that the total distribution was an addition of 



two separate component Pareto laws, called a "Bi-Pa- 
reto law," and expressed by its probability density func- 
tion as follows: 



fir) — p + a — —■ — 



(1) 



where p and q are normalization constants for the sum- 
mation of both Pareto component functions to integrate 
to unity and are such that/? + q = 1, so that (1) actually 
depends on five parameters (a, b, a, p, and either p or 
q). 

Note that to perform a fit to the histograms of the 
random variable log 10 (r) and not t itself, as shown in 
Figs. 2 and 3, the r distribution must first be multiplied 
by r. Also note that the fit criterion used throughout 
this study was to perform a least squares fit — that is, 
minimize the standard deviation between the logarithms 
of the measured and "theoretical" distributions. 

The standard deviation was found to be 0.13 1 for this 
Bi-Pareto model. The first of the two components of the 
model represents the shorter time intervals and what 
Jaappens within a rainy episode. As mentioned earlier, 
this component will be strongly dependent, through its 
a and a parameters, on the size of the measuring cross 
section of the instrument. The second component, which 
accounts for the larger time intervals, expresses the 
spacing between rainy episodes and will not depend on 
the instrument size. 

The least squares fit of this function to the histogram 
of Fig. 3a after applying a normalization constant C 
[since, for simplicity, it was chosen to express histo- 
grams as "number of events per bin," that is, the number 
of log 10 (T) samples found between two consecutive 
thresholds] gives the curve plotted as a solid line, with 
the following model parameters for r in milliseconds: 



[a = 133, 
[a = 1.35, 



b = 1220, 
/3 = 0.68, 



q = 1.12 X lO" 2 , 
P = 1 ~ o, 



= 1.235 X 10-«J. 



(2) 



From the point of view of goodness of fit to exper- 
imental data, this model may be considered excellent. 
However, other mathematical expressions could as well 
be used instead of (1) with a similar quality of fit. Sim- 
plicity and the possibility to separate short-term from 
long-term events have decided in favor of model (1). 
However, this model nevertheless raises serious prob- 
lems. The mean value of the interdrop time interval, 
provided it exists, is given by 



Jo 



aa a r 
(T+a)° +i 



dr + q — dr. 

q Jo (T + *)*♦' 

(3) 

It exists only if a or /3 are both greater than unity. 
However, here, while a is greater than 1, this is not the 
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Fig. 4. Distributions of interdrop time intervals obtained from 4 
months of data (O) and from 25 months of data (+). 



case for j8. Therefore, applying the model would imply 
that 

• the process is not stationary and that 

• the process has an "end," meaning that the average 
number of drops N over an interval [t, t + 8] as- 
ymptotically decreases as follows: 



E[N(t. t + 5)] -> 0 
for / oo 



V 8 > 0, 



which amounts to predicting an endless and overwhelm- 
ing drought at some time in the future (Feller 1957; 
Karlin and Taylor 1975). Obviously, none of these con- 
sequences are physically acceptable. To remove these 
inconsistencies, it is necessary to admit that the prob- 
ability density function f{r) of (1) becomes invalid for 
very large r's. 

To investigate this point, the low-resolution dataset, 
which covered a full noninterrupted 25-month period, 
was used. Although this dataset has a lower time res- 
olution, the very large time period it covers allows lon- 
ger interdrop time intervals to be studied. It leads to a 
histogram of its own, which gives an indication of how 
the tail of the overall distribution should look. We chose 
to display it on the same graph as that obtained from 
the high time resolution dataset. To combine the short- 
and long-term histograms, the latter was normalized in 
such a way that it appears as a continuation of the former 
at the tail of the high-resolution dataset since what we 
are looking for is essentially the behavior of the his- 
togram at longer interdrop time intervals; that is, we are 
checking that the histogram actually falls down for large 
t's. The result is shown in Fig. 4, in which the circles 
and crosses identify the short- and long-term data, re- 
spectively. This pseudonormalization procedure leads to 
counts less than unity in the second histogram for ab- 
scissa greater than about 3 days. What is important to 



note in this figure is an abrupt decrease in the histogram, 
somewhere around a 1-day interdrop interval. The ex- 
perimental combined histograms thus obtained suggest 
that an exponential decrease occurs for the larger in- 
terdrop time intervals. 

To take this exponential decrease into account in the 
model, it is, however, not possible to directly multiply 
the Bi-Pareto distribution function f{r) by an exponen- 
tially decreasing factor of the form e~ x ° T , as it would 
not integrate to unity. This problem can be overcome 
when noting that 



e-*orf ( T ) dr = 1 - 
Jo Jo 



A 0 F c (T)e-V dr, 



where F c (t) is the complementary cumulative distri- 
bution function defined as F c (r) = 1 — \ T 0 f(x) dx. From 
this, it may be seen that to have a statistical distribution 
whose main form is f(r) and that shows an exponential 
decrease for large t's, the Bi-Pareto distribution f(r) 
has to be transformed as follows: 

<p( T ) = e -^[f{r) + A 0 F c (r)]. (4) 

For sufficiently small values of A 0 and for small r's, 
<P(t) 888 /(t), while for large r's, <p(r) = e~ A ° T /(r), which 
shows that the function has the desired exponential de- 
crease. When applied to the Bi-Pareto law (1), the trans- 
formation suggested in (4) leads to the following prob- 
ability density function: 



f(r) = e~ A ° T 



a"[a + A 0 (r + a)] 
(r + a) a+l 

b*[p + A 0 (t + b)] 



+ q 



(t ■+ by +l 



(5) 



which will be called the Complete Bi-Pareto model, and 
to the associated cumulative distribution function 



F(t) = 1 - e~ x * T 



a a 



b? 



(t + a) a 72 (t + by 



. (6) 



Keeping the values of the parameters already obtained 
above, the value of A 0 can be derived through a best fit 
to the experimental histogram of Fig. 4, again using a 
least squares fit between the logarithms of the model in 
(5) (after multiplying it by t, as mentioned above) and 
the experimental histogram. The value found was A 0 = 
10~ 8 . The theoretical curve thus obtained is plotted as 
a solid line in Fig. 4, where it may be seen that the 
probability density function (5) (after multiplication by 
an appropriate normalization constant) is in very good 
agreement with the data collected during 25 months. 

This distribution, as expected, provides a finite av- 
erage interdrop time .interval. Since a\ 0 and b\ 0 are 
small, the average interdrop interval fi (3) may be given 
a simple asymptotic expression 



fJL = <T> * p 



a - 1 



+ gb^rm - p), (7) 
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Fig. 5. Occurrence rate function s(r). 



which leads to the numerical result that fx = 0.38 + 
1 A = 1.78s for our experimental setup and observation 
period. The value of the occurrence, ra.to, ajs.0, known ^ 
the hazard function , s(t) = '/(r)/[l - F(t)], which in- 
dicat es the probability tW a WiHdrop, arrived be t w een 
time points r and r + dr, provided none has arrived 
before time r, is shown in Fig. 5. For small t's, s(t) is 
constant. From 0.1 s to 1 day it decreases as T/t, which 
may be interpreted by saying that the more one has been 
waiting, the more one will be waiting. In other words, 
the more the rain rate decreases, the more it is likely to 
end. This is a common sense observation that leads to 
the subjective feeling of clustering. Fortunately enough, 
this trend does not last, as s(t) asymptotically goes to 
A 0 , which is a physical expression of the revisited prov- 
erb: "After sunshine comes the rain!" 



b. Relationship between size and time of arrival 
distributions 

The question of how the survival law defined in (5) 
is generated will now be addressed. One of the re- 
markable properties of Pareto distributions is that they 
can be obtained by a combination of an exponential 
distribution and a gamma distribution (Harris 1968). 
Assuming these two distributions write 



p(r) = Xe~ 



(8) 



and 



p(A; a, a, A 0 ) 



T(a) 
0 



(A - Ao)- 1 *-^, A > A 0 , 



A < A 0 



(9) 



their mixture leads to the Complete Bi-Pareto's distri- 
bution, namely, one has 



f 



■II'(A)p(t) d\ 



fl«[q»A 0 (T ■+.<!)] , A 0 (T + 6)] 



(r + a) a 



(t + by +i 



(10) 

which is the expression of the Complete Bi-Pareto mod- 
el of (5). In this expression, 11(A) is a combination of 
two "shifted" gamma distributions as defined in (9): 

n(A) = pp(A; a, a, A 0 ) + gp(\; ft b, A 0 ). (11) 

This raises the question of what such a distribution 
might physically correspond to. First looking at distri- 
bution (8), A represents the inverse mean time between 
two raindrops or the mean number of drops per time 
unit. Over a time period T, AT then asymptotically rep- 
resents the number of drops with "parameter" A. Since 
the drop size is the main characterizing parameter, it is 
therefore tempting to link A with the effective drop di- 
ameter D, the involved relation obviously being mo- 
notonous. 

This assumption can be tested using many different 
distributions that relate A with D. It is desirable to pre- 
vent the drop size distribution model from overesti- 
mating very small raindrops, which are the most nu- 
merous, as this might strongly bias physical quantities 
derived from the model, while avoiding singular points. 
After several tests, the following \-D relation was 
adopted and will be justified below: 

A = A - A 0 - K(D - D,)"[l - e-^-oo"]. (12) 

The following drop size probability density function is 
then obtained: 



Yl(D) = 11(A) 



dk 



dD 



be 



T(a) 



vA 



— + vkD v - x {KD v - A) 



(13) 



Since (13) is a probability density function, it can be 
compared with the experimental histogram (Fig. 3b) 
only after having been normalized to the total number 
of drops. For experimental reasons, drops with diame- 
ters under 0.3 mm were rejected from the dataset as 
dubious, so that parameter £>, in relation (13) is, of 
course, impossible to determine through a direct fitting 
process to the available data. 

When £>, is set to different values below 0.3 mm, the 
parameters indicated in Table 1 are found for the drop 
size distribution model. The influence of the model on 
such a physical quantity as the mean drop volume is 
indicated in the (v) column of Table 1. 
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Table I. Best-fit parameters in relation (13) (C is the 
normalization constant found in the fitting process). 













» 






(mm) 


V 


K 


k 


a 


(mm 3 ) 


C 




0.05 


0.55 


3.5 X 10- 2 


1.66 


0.062 


0,14 


4.05 X 


10 s 


0.1 


0.56 


3.4 X 10" 2 


1.70 


0.061 


0.18 


3.65 X 


1 0 s 


0.2 


0.57 


3.3 X 10- 2 


1.70 


0.060 


0.25 


3.04 X 


10 s 



The obtained fitting results are all very good (standard 
deviations a close to 0.06 count units), as shown in Fig. 
6, so that this criterion alone cannot be reliably resorted 
to here. This result is rather interesting for two reasons. 

1) The final model that has been arrived at provides 
a coherent view of the survival and drop size distri- 
butions. Nevertheless, this does not mean that a "rea- 
son" why these distributions have the shape they have 
has been found. However, the fact that the model implies 
a coupling between drop sizes and arrival times opens 
new perspectives. 

2) The value found for v is the most surprising result. 
When referring to the work of Atlas and Ulrich (1977), 
who suggest the relation between the fall velocity at 
ground level and diameter 

V(D) = 3.78D 0 - 67 , (14) 

where D is expressed in millimeters, it may be noted 
that parameter v that has been found and is indicated 
in Table 1 differs only by 15% from the value of 0.67 
shown in (14). The large value of k leads to the result 
that whenever D is greater than 0.5 mm, (12) can be 
approximated by K(D - D x ) v and is thus found to be 
very close to relation (14). The effect of the 1 - 
e -Ko-Do v term j s mat n(Z>) approaches 0 when D ap- 
proaches D,, which is satisfactory. Here again we can 
regret only the lack of reliable data for small raindrops. 
Therefore, everything tends to show that the quantity 
A in formula (12) is actually proportional to the ground- 
level fall velocity of raindrops. 

c. Applications of the proposed model 

Rainfall rate is universally employed to characterize 
rain. It would therefore be interesting to investigate what 
the present model would be able to reveal in this respect. 
This question does not find a straightforward answer, 
as it is related to the counting aspect of the process, 
which hardly lends itself to analytical computations, at 
least in the case of survival law (5). However, classical 
calculations (appendix B) lead to a dual Laplace trans- 
form of the probability that a precipitated water height 
h be exceeded over a given [0, /] time interval 



QM s) = 



1 - f(s) 1 - p(u) 



JJLS 2 U 



1 - p(u)f(sY 



(15) 



where p(u) is given by (B9) and p, is the average in- 
terdrop time given by (7). Expression (15) allows the 
two following problems to be tackled easily. They are 




Fig. 6. Experimental (October 1993 to February 1994) and theoret- 
ical drop size distributions. 



selected not only to check the model's validity but also 
to show some of its practical applications. 

1) Average precipitated water height 

The average precipitated water height can be obtained 
either directly from the basic renewal theorem applied 
to a cumulative process (Feller 1957) or through the 
relationship lim u _> 0 Q c (u, t) = Jo Q c (h, t) dh, which leads 
to 



WO) = ft, 



(16) 



where (v) is the average raindrop volume and A is the 
collecting area of the instrument. When using the value 
of ji reported above — that is, ji - 1.785 s — and size 
distribution (13), which leads to (v) = 0.26 mm 3 , the 
average yearly precipitated water height is 

<#(0> y e a r = 549 mm. 

This value is in good agreement with that measured by 
means of a standard rain gauge located on the same site 
as the disdrometer, which led to a total rain height of 
563 mm for the period from May 1992 to April 1993. 
However, this agreement should be regarded with cau- 
tion, as the variability of cumulative rain height from 
one year to the next may exceed 100 mm. 

Can this model explain such a variability? Appendix 
C gives a computation of the main elements that are 
needed to estimate the standard deviation of rain height 
that should be expected. Numerically, it is found to be 
o- H l{H) = 3%. Such a small value is indicative of the 
fact that the model requires further refinement to account 
for seasonal or climatic variability. This point will be 
addressed in a further paper when sufficient data are 
available. 
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Fig. 7. Comparison between asymptotic and best-fit forms of the 
Laplace transform f{s) of /(t), as defined by (19) and (20), respec- 
tively. 



2) Conversion between rain rates obtained 
with various integration time 

Designers of radio systems use predictions based on 
rainfall statistics derived from rain gauges with 1-min 
integration times (ITU-R 1990). However, the integra- 
tion times of rain gauges may differ from one place to 
another and may reach 10 min or even an hour. In such 
cases, it is useful to derive statistics that would have 
been obtained for 1-min integration times from those 
available at other integration times. Traditionally (Segal 
1986), one uses a conversion factor that depends on the 
considered exceedance probability and on the rain- 
gauge integration time. 

The present model allows this problem to be ad- 
dressed in a different way. When focusing on higher 
rain rates, the asymptotic value of p(u) for u -» 0 may 
be used in Q c (u, s); that is, p(u) = 1 — ({v)/A)u = 1 
- Bu. This leads to (cf. appendix B) 



i 



u + 



i - m 

Bf{s) 



, (17) 



with the following inverse Laplace transform in u: 

Q£h, s) ~ 1 - ^ e -M[i-/(#[W]) (18) 

fJLS 2 f(s) 

For integration times greater than 1 min, the inverse 
transform in (18) will be influenced only by values of 
s that are greater than 10~ 4 . Thus, the asymptotic form 
of f(s) may be used: 



for s ranging from 10~ 8 to 10" 3 , by a function of the 
following kind: 

f(s) « 1 - ds*, (20) 

where £ = 0.92 and d = 10 233 . Substituting (20) into 
(18), the following equation is obtained: 

QXK s) = - s ^e-^\ (21) 
/* 

Applying the scaling property of the Laplace transform 
leads to 



QM kt) 
from which it is derived that 



K4 



e \h' ') = kt ~ lQAh - kt) > 



As) - 1 - P 



as 
a- 1 



qb*s{s + A 0 )^T(1 - 0). (19) 



or using rain rate r measured over interval [0, t], with 
r = hit, 

Qdr, t) = k^Q c {rki~\ktl (22) 

where g c (r, t) is nothing but the cumulative probability 
function of rain rate r that would be obtained with a 
rain gauge having an integration time t. Formula (22) 
thus allows a rain-rate distribution observed using an 
integration time kt to be changed into one that would 
have prevailed for an integration time t. 

To show how this model can be applied, data obtained 
using the same disdrometer, during a 1-yr period not 
overlapping the one used to develop the interdrop time 
distribution model (5), were analyzed using 1- and 10- 
min integration times. 

The ITU-R model mentioned above has been derived 
from measured cumulative distributions, for the 5- and 
10-min to 1-min conversions, and cannot be extended 
to other integration times. It has been obtained by fitting 
the ratio of the 5-min to the 1-min cumulative distri- 
butions measured in Canada, without considering the 
fine structure of the rain's temporal distribution. 

This temporal information is contained in the present 
Complete Bi-Pareto model, which has led to the inte- 
gration time conversion model suggested in (22). To 
check the validity of the present model and compare it 
to the ITU one, it is interesting to apply it to the common 
10-min to 1-min conversion problem. 

Under these conditions, factor k in (22) is equal to 
10, and k*~ l is equal to 0.768, assuming £ = 0.885, as 
suggested in the previous section. From this, if P, is the 
probability that a rain rate r,, measured with a 1-min 
integration time, is exceeded, then the probability P, 0 
that a rain rate r 10 , measured with a 10-min integration 
time, is exceeded, would be given by 

PM = ^(r.VO.768 

and 



Figure 7 shows that this function may be approximated, 



r 10 = r, X 0.768. 
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Fig. 8. Comparison between measured cumulative distributions of 
rain rate with 1- and 10-min integration times, given by the present 
model [(22)] and the ITU-R model. 



In the ITU-R model, the corresponding relation is as 
follows: 

= P,(r,) 

and 

r x X 0.836 
ri ° ~ [Px(r x ) X 100] 00736 ' 

The 0.836 factor in the above expression was defined 
as corresponding to climatic zone E, the site where our 
measurements were performed. 

Figure 8 shows the cumulative distribution of rain 
rate using integration times of 1 min (continuous line) 
and 10-min (circles). The distribution obtained accord- 
ing to the present model [Eq. (22)] is plotted (dash-dot 
curve), and a comparison is provided with the ITU-R 
model (dots). Not only does this figure show a very 
good agreement between our model and observations 
but the present result is also in good agreement, at least 
up the 50 mm h _1 level, with the ITU-R model obtained 
using a quite different and purely empirical approach 
and based on a very different data corpus. 

The advantage of this model is that it is universal, in 
the sense that any combination of integration times can 
be used — contrary to the ITU one — that imposes 5-, 
10-, and 15-min integration times. 

5. A fractal point of view 

Atmospheric phenomena, like many natural phenom- 
ena, have self-similarity properties in that their form 
does not change whatever the scale at which they are 
observed. Lovejoy and Mandelbrot (1985) presented a 
rain model with only one fractal dimension (which is 
the parameter relating the statistical properties to all 
observation scales). However, this approach does not 



allow all observations to be represented. Moreover, the 
fractal dimension seems to depend on the observation 
techniques used (e.g., rain gauge, radar, etc.). In a recent 
paper, Olsson et al. (1993) analyzed the similarity prop- 
erties of a time series of data from a very high resolution 
rain gauge (0.035 mm min -1 ). It is interesting to look 
at what our data and model may lead to using the same 
approach. The box-counting method developed by these 
authors was chosen due to its greater simplicity than the 
other suggested analysis methods. 

The total observation time T is divided into n con- 
tiguous intervals of length L. Thus, T = nL. The total 
number of occupied intervals N(L) y that is, those in 
which at least one raindrop has been observed, is then 
counted. If the data form a one-dimension fractal, then 
we must have 

N(L) a L-y, (23) 

where y is defined as the fractal dimension of the pro- 
cess. When plotting N(L) as a function of log(L), a 
straight line with slope —7 is obtained. Of course, when 
L->T, N(L) -» 1, and, asymptotically, one should have 

Mm N(L) = n = TIL. 

L->T 

Similarly, if J, is smaller than the minimum time sep- 
ara tion between raindrops. N(M = m. where m is the 
toTai number or' drops observed. 

In any case, the behavior represented by (23) can 
occur only far from the two ends of the total observation 
period. Figure 9a shows the result of such an analysis 
when applied to the longest available period (4 months) 
of high-resolution data. 

Except for the asymptotic regions mentioned above, 
the curve found (circles in Fig. 9a) is very close to a 
straight line with a slope, — y, which has been found to 
be -0.82 and to prevail in the range from about 10~ 2 
to 10 4 min. When comparing these results to those found 
by Olsson et al. (1993), from which Fig. 9b was taken, 
a fundamental difference may be noted. Namely, in our 
case, a single fractal dimension is observed up to 10 4 
min, whereas the use of a 1-min time resolution tipping- 
bucket rain gauge leads to a three-segment curve with 
slopes of -0.82, -0.37, and -1 for increasing box 
lengths L. A first change in slope occurs near the 45- 
min abscissa. 

The equality of the obtained fractal dimensions 
(—0.82) for smaller L suggests that the second segment 
obtained by Olsson et al. may be an instrumental effect, 
as the integration inherent to the measurement principle 
of rain gauges has a strong influence on the obtained 
data. This point of view is corroborated by the influence 
of the threshold value on the fractal dimension (Olsson 
et al. 1993). 

Assuming the Complete Bi-Pareto model provides an 
accurate representation of the temporal behavior of rain, 
it should then be possible to directly infer the curve in 
Fig. 9a from the model. Thus, let an arbitrary time in- 
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Fig. 9. (a) Box-counting distributions N(L) for the 4-month data 

(O) and derived from the Complete Bi-Pareto model ( ). (b) 

Experimental box-counting distribution N(L) for 2-yr measurements 
(O) with best fit ( ) from Olsson et al. (1993). 



terval be of length L. Then the probability that it is 
"occupied" is equal to the probability that at least one 
drop falls into this interval, or, in other words, is the 
value of the cumulative distribution function of the re- 
sidual lifetime taken at time L. Using the cumulative 
distribution (6), this leads to 



-f 

Jo 

■f 

Jo 



1 - F{r) 



dr 



£-A 0 T 


pa a 




(r + a) a 



qb* 



(t + by 



dr. 



If the total observation time Tis divided into n intervals 
of length L, the theory of large numbers may be applied 
for estimating the mean value of the occupied intervals: 



(N(L)) = -G(I), 



or 



(N(L)) 



qb* 



(t + a) a (t + by 



dr. 
(24) 

Expression (24) is either evaluated by means of in- 
complete gamma functions or directly through a simple 
quadrature. The result obtained based on the parameters 
(2) found for the Complete Bi-Pareto model is plotted 
in Fig. 9a (solid curve). 

A very good agreement can be noticed between the 
theoretical and experimental curves except for a slight 
vertical offset. This offset is likely due to an insuffi- 
ciently accurate estimation of the average interdrop in- 
terval fx, which in turn is derived from parameter A 0 
that has been shown to be known only approximately. 
This result corroborates the assumption of independence 
of raindrop times of arrival and therefore strengthens 
the validity of our modeling in terms of renewal pro- 
cesses. 

It should also be underlined that the proposed model 
is richer than the fractal approach alone, which may also 
be suspected to be very sensitive to thresholding. 

6. Conclusions 

The use of a high-resolution disdrometer allowed a 
raindrop time-of-arrival time series to be obtained in 
combination with their sizes. An analysis of this time 
series over a total duration of 14 months, with a max- 
imum of 4 continuous months, led to the following con- 
clusions. 

• The time series can be interpreted as the result of a 
renewal process whose survival law is of the Com- 
plete Bi-Pareto type. 

• The survival law is related to the size distribution. 

The interest of the model, in addition to a very good 
agreement with the collected data, also resides in its 
ability to explain various aspects (e.g., precipitation in- 
tensity and fractal behavior) of rain. 

Two applications useful in the fields of meteorology 
and telecommunications have been shown and provide 
a good validation of the model. 

Some particular aspects however still deserve atten- 
tion. 

• The influence of small raindrops currently incorrectly 
measured down to a minimum diameter of 0.3 mm. 
Although such an accuracy may at first seem to be an 
achievement for this kind of apparatus, the measured 
distribution curves and derived models show that mea- 
suring even smaller diameters should be made pos- 
sible. 

• The model's universality. Can rain be described by 
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Fig. 10. Effect of diameter selection threshold on interdrop time 
distribution (3-week period). 



means of the same model anywhere in the world? Of 
course, nothing allows such a conclusion to be drawn. 
However, the fact that the same fractal dimension is 
found in Paris and Sweden is encouraging. 

Other more fundamental aspects, such as the com- 
bination of spatial and temporal effects, were not ana- 
lyzed here and would allow further applications of the 
model to be envisioned in various fields, such as tele- 
communications, etc. 

Nevertheless, the suggested model seems to be of 
sufficient interest that further experiments, it is hoped, 
may be carried out in other climatic conditions and over 
longer uninterrupted time periods in order to find dif- 
ferent values for its climate-dependent parameters. 
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and therefore different time period than those that led 
to a, b, a, fi, and q. Therefore, its value should be 
considered only as an indication. Nevertheless, this 
model appropriately describes observations and this, 
over long durations, was far from obvious at first sight. 

However, it may be wondered whether rejecting rain- 
drops with sizes less than 0.3 mm from interdrop in- 
terval data as well as from drop diameter data, as was 
done here, had any effect on the model's shape. The 
model itself provides an answer to this question, as one 
of its key properties is a mixture between a size (or 
velocity) distribution and an exponential survival law. 
Accordingly, the renewal process can be divided into 
subprocesses corresponding to different drop size class- 
es. A classical computation (e.g., Feller 1957) shows 
that if drops are selected with sizes that, for example, 
lie within a range / = [D mini D^], a new regenerative 
process is generated, the distribution function of which 
has the following Laplace transform over an interval /: 



m = 



f(s)u 



1 + (u - !)/(*)' 



where u = Prob[D min < D < D ma J. 

For times ranging from a few minutes to a few hours, 
f(s) can be read as f(s) = 1 — g(s) with g(s) 1 and 
Si s )s->o ~~* u - Under the same conditions, the following 
result may be obtained: 

m ~ 1 - 8 -f- 

Since the shape of /(r), for r -» °°, is directly related 
to that of its Laplace transform, for s -» 0, it may be 
expected that the slope of //(r), in logarithmic coor- 
dinates, for large r's, remains unchanged whatever the 
chosen / interval. Figure 10, which gives the histograms 
corresponding to four drop size selection thresholds [D 
> 0.5 mm (O), D > 1 mm (+) and D > 2 mm (X)], 
is in very good agreement with this conclusion except 
for the fluctuations in the tail portion of the distributions, 
which are due to the reduced number of drops after 
selection. 

It therefore appears that the model was not changed 
by eliminating drops smaller than 0.3 mm from the da- 
taset. 



appendix A 

Effect of a Diameter Threshold for Selecting 
Analyzed Data 

In section 4a we showed how to derive a joint rep- 
resentation of experimental distributions of interarrival 
times and drop sizes. The only remaining "free" pa- 
rameter was Z),. For the time being, it may not be de- 
rived, as it is apparently related to the disdrometer's 
sensitivity level. In addition, parameter A 0 was estimated 
only from long-term data that covered a much longer 



APPENDIX B 

Computation of the Precipitated Water Height 

Let us consider the time interval T„ between n rain- 
drops and determine its distribution, with an arbitrary 
time origin: 

r„=T r + §T ( . (Bi) 

The t/s are independent, identically distributed ran- 
dom variables with a probability distribution f(r) given 



August 1998 



LAVERGNAT AND GOLE 



817 



by (7). The variable r r is also independent of t, and is 
the residual survival time whose probability law is given 
by 



A(Tr) = 



1 ~ F(T r ) 



(B2) 



Thus, T n follows a probability density function that 
is obtained through multiple convolution, namely, 



Kit) = /,*/<-><*> 



(B3) 



Taking its Laplace transform (keeping the same no- 
tations for functions and their Laplace transforms), 
where s is the variable in the Laplace space, and using 
obvious notations leads to 



K(s) = mm- 1 = - — — 

jJLS 



(B4) 



If R(t) is the number of drops fallen during interval 
[0, t] 9 the probability that it is less than n is 

?r{R(t) < n) = ?v(T n > t) = 1 - K„(t), (B5) 

where K n is the cumulative distribution function cor- 
responding to density function (B3): namely, 



KM = \ K{x) dx 
Jo 



and, of course, K 0 (f) = 1. From (B5), the probability 
U n (t) that R(t) = n can be defined as 

11.(0 = Pr{J?(0 = n} 

= ?v{R(t) > n) - Pr{/?(0 > n - 1} 

= K n (t) - K n _ x {t\ (B6) 

where R(t) is an integer random variable with a gen- 
erator function g(z, t) given by 

00 

g(z, t) = 2 n«z". (B7) 
Laplace transforming (B7) over variable t leads to 

Since /^(s) = £„(.s)/s and K 0 (s) = 1/s, 



g(z. 5) = - + 

5 5 



which, after summation of the geometrical series, leads 
to 



s s 1 - zf(s) 



(B8) 



Besides, a raindrop of diameter D and volume v = 7rZ) 3 / 
6, falling through a collecting surface area A, contributes 
a water height h = v/A = wD 3 /6A. The probability 
density function of h is derived from that of D since 



P<*> = P(D) § 



Its Laplace transform is directly obtained as 

J "00 Too 
p(h)e- h "dh e~W^p(P)dD. (B9) 
0 Jo 



Actually, the total water height //(/), fallen over in- 
terval [0, f], is a cumulative process since in our model 
the drop sizes are independent of arrival times. Using 
the total probabilities theorem, the probability density 
function of H can then be written 

q(h, t) = ?r{h <H^h + dh) 

= 2 n„(r)Pr{A < H < h + d&U^.). 
Applying a Laplace transform over variable h gives 

9(% o = 2 n B (f)p-(«). (bio) 

A comparison of (BIO) and (B7) immediately shows 
that a second Laplace transform over variable t gives 
q(u, s) = g[p(w), s] or, explicitly, 

1 /,(!) 1 ~ p{u) 

q{u, s) = - + — (Bll) 

5 5 1 - p(u)f(s) 

A dual Laplace transform of the cumulative distribution 
function can then be written as 

0( 1 Ad) 1 ~ P(u) 

' su su 1 - p(u)f(sY 

The complementary distribution Q c (h, t) = Pr{H > h, 
t} has a dual Laplace transform of the form 

QM s)= l ~ & - 1 - ■ (B12) 

1 - p(u)f(s) 

Theoretically, the problem is "solved." Unfortunate- 
ly, it was not possible to find an analytical inversion for 
(B12). 

APPENDIX C 

Variability of the Precipitated Water Height 

Let./V be the number of drops counted_aver_inter-val 
[0, /]. Each drop / contributes a water height h iy which 
is independent and identically distributed with a distri- 
bution p(h). The total water height fallen during the 
same time interval is given by 



H=2 K 



(CI) 



Thus, H appears to be a random sum of independent 
and identically distributed random variables. A conven- 
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tional computation (e.g., using Laplace transforms) 
gives the average value, which is nothing but (16), and 
the variance a 2 = [(N 2 ) - (N) 2 ](h) 2 + [(h 2 ) - 
(h) 2 ](N\ or 

crl= <r 2 M 2 .+ <rl(N). (C2) 

Here, (h) 2 and v 2 h can be computed easily from the drop 
size distribution, and (N) = tf/x and a N can be appro- 
priately approximated by using (Karlin and Taylor 1 975) 
a 2 N ~ tcr 2 /^, where a 2 is the survival law's standard 
deviation. A computation involving a Laplace transform 
shows that for our model 

(t 2 ) - 2/?f + 2qb^\ - /3)Ag- 2 (C3) 

or, with the numerical values found, (r 2 ) = 6.1 X 10 10 . 
Thus, a — 2.5 X 10 5 , the second term of (C2), is neg- 
ligibly small relative to the first one, and the standard 
deviation of H may be written 

(2\ ,/2 
W • (C4) 

During a year, (C4) leads to <j h I{H) 3%, and during 
a period of 1 month to 
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